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HIGHLIGHTS 


•  Two  mathematical  models  of  lithium-ion  cell  with  blended-electrodes  are  developed. 

•  Equilibrium  model  can  easily  calculate  equilibrium  characteristics  of  a  blended-electrode. 

•  Dynamic  model  can  predict  dynamic  performance  of  a  lithium-ion  cell  with  blended-electrodes. 

•  Developed  model  is  validated  with  experiment  data  and  shows  good  agreements. 
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Many  lithium-ion  batteries  have  electrodes  made  of  multiple  types  of  active  materials  in  order  to  take 
advantage  of  each  active  material.  According  to  the  types  and  blending  ratios  of  constituting  active 
materials,  there  can  be  numerous  combinations  of  blended-electrodes,  resulting  in  different  performance 
characteristics.  This  paper  introduces  how  to  numerically  model  lithium-ion  batteries  with  blended- 
electrodes,  and  predict  their  performance.  Further,  this  paper  demonstrates  how  to  apply  the  devel¬ 
oped  model  to  designing  new  batteries.  First,  we  propose  an  equilibrium  model  which  can  predict 
equilibrium  characteristics  of  a  blended-electrode.  Second,  a  physics-based  lithium-ion  cell  model  with 
blended-electrodes  is  proposed.  Developed  model  is  validated  with  experimental  data  and  shows  good 
agreements. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Although  lithium-ion  batteries  commercialized  by  SONY  have 
several  merits  such  as  larger  capacity,  lighter  weight,  and  good 
coulombic  efficiency  compared  to  previously  used  batteries  of 
other  chemistry  such  as  nickel-metal  hydride  batteries,  there  has 
been  doubt  that  lithium-ion  battery  is  not  safe  enough  to  be  used 
for  automotive  power  source.  However,  recently,  lithium-ion  bat¬ 
tery  is  getting  much  more  public,  industrial  attention.  Several  au¬ 
tomakers  have  launched  hybrid  electric  vehicles  (HEV),  plug-in 
hybrid  electric  vehicles  (PHEV),  and  battery  electric  vehicles  (BEV) 
powered  by  lithium-ion  batteries  after  GM  successfully  launched 
the  first  commercialized  plug-in  hybrid  vehicle,  VOLT™  which 
loads  LG  Chem’s  high  performance  lithium-ion  batteries.  In 
contrast  to  small  batteries  for  portable  electronics,  automotive 
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lithium-ion  batteries  must  sustain  severe  operating  conditions  and 
requirements.  Therefore,  it  requires  much  more  time  and  efforts  to 
develop  large-formatted  automotive  batteries  than  small  batteries 
for  mobile  applications. 

Typical  lithium-ion  batteries  are  composed  of  negative  active 
materials,  positive  active  materials,  electrolyte,  separator,  additives, 
and  current  collectors.  Active  materials  are  the  most  important 
among  them  since  they  determine  power  and  energy  of  a  battery. 
Many  negative  and  positive  active  materials  have  been  searched  up 
to  now  and  each  of  active  material  has  its  own  characteristics  such 
as  operating  potential  profile,  capacity  and  power.  For  example, 
LiMn204  has  good  rate  capability  and  high  power  with  relatively 
low  cost  and  less  safety  problem.  But,  it  shows  low  capacity  and 
Manganese  (Mn)  dissolution  problem  also.  In  contrast,  Li(Ni,  Mn, 
Co)02  has  more  capacity  compared  to  LiMn204,  but  shows  lower 
power  due  to  the  layered  structure  and  lower  operating  potential 
window.  Olivine  material  such  as  LiFeP04  is  considered  as  one  of 
promising  positive  materials  due  to  low  cost,  improved  safety  with 
large  capacity.  Many  researches  to  overcome  low  electronic 
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conductivity  and  very  low  lithium  diffusion  rate  of  LiFeP04  (briefly 
‘LFP’)  particle  have  been  conducted  [1-4].  As  briefly  listed  above, 
most  of  active  material  have  both  merits  and  demerits  together. 
Therefore,  it  is  important  to  select  proper  active  materials  for 
specified  purposes.  Since  large  battery  capacity  is  the  most 
important  factor  for  full  electric  vehicles  (EV),  Li(Nii/3Mni/3Coi/3)02 
(briefly  ‘NMC)  can  be  selected  as  positive  active  material  and 
graphite  can  be  used  for  negative  electrode.  In  contrast,  high  power 
is  the  most  required  quality  for  hybrid  electric  cars  (FIEV).  In  this 
case,  spinel  materials  such  as  LiMn204  (briefly  ‘LMO’)  and  amor¬ 
phous  materials  such  as  carbon  can  be  adopted  for  positive  elec¬ 
trode  and  negative  electrode,  respectively.  When  both  capacity  and 
power  are  important,  multiple  active  materials  can  be  blended 
together,  so  called  ‘blended-electrode’. 

Several  papers  about  experimental  work  on  blended-electrodes 
have  been  published.  Mitigating  Manganese  dissolution  and 
improving  capacity  retention  by  blending  different  active  mate¬ 
rials  together  were  considered  important  by  several  researchers 
[5-7].  Fergus  [8  extensively  reviewed  cathode  materials  including 
composite  electrodes.  He  categorized  electrode  with  different 
mixed  particles  and  electrodes  with  coated  particles  layered  as 
composite  electrode.  Whitacre  et  al.  [9]  constructed  fully  mixed 
(blended)  electrodes,  fully  segregated  electrodes  and  layered 
electrodes  by  using  LiFeP04  and  Li(Li0.i7Mn0.58Nio.25)02.  They  tried 
to  take  advantages  of  composite  electrode  from  each  active  ma¬ 
terial  (high  rate  capability  of  LiFeP04  and  large  capacity  of  Li(L- 
io.i7Mno.58Nio.25)C>2  and  they  found  that  segregated  active  material 
configuration  is  the  most  promising  solution  with  those  two  active 
materials. 

As  briefly  reviewed  above,  it  is  worthwhile  to  consider  blending 
several  active  materials  together  to  improve  cell  performance. 
However,  blending  different  active  materials  may  cause  some 
problems  in  a  complex  slurry  mixing  process  or  electrode  coating 
process,  which  may  result  in  increased  manufacturing  cost.  In 
addition,  much  time  and  cost  would  be  required  to  find  optimum 
blending  ratio  for  target  performance  of  a  cell.  Actually,  it  takes 
much  time  and  cost  to  finalize  cell  chemistry  because  cell  makers 
go  through  several  screening  tests  for  candidate  materials.  Further, 
numerous  new  active  materials  are  being  introduced  by  material 
makers.  Although,  it  is  good  for  cell  makers  to  have  many  options  to 
choose,  required  time  and  cost  for  a  new  cell  development  may 
increase  as  the  combination  number  of  materials  grows  larger.  If 
the  performance  of  a  new  cell  with  blended-electrodes  can  be 
mathematically  predicted  before  the  actual  cell  is  constructed,  cell 
development  cost  will  be  significantly  saved. 

Few  publications  about  modeling  work  on  blended  electrodes 
can  be  found  up  to  now.  Darling  and  Newman  [10]  started  a 
modeling  work  on  the  electrode  with  two  distinct  particle  sizes. 
Gomadam  et  al.  [11  ]  developed  a  mathematical  model  simulating  a 
composite  electrode  made  of  carbon  monofluoride  (CFX)  and  silver 
vanadium  oxide  (SVO)  for  medical  applications.  Although  their 
model  is  relatively  simplified  with  many  assumptions,  their 
modeling  result  showed  good  agreement  with  the  experimental 
result  from  a  small  half  cell  test  within  moderate  C-rate  discharging 
condition.  Albertus  et  al.  [15-17]  conducted  more  extensive 
modeling  work  on  blended  electrodes.  They  modified  their  math¬ 
ematical  cell  model  to  treat  electrodes  made  of  multiple  types  of 
active  particles,  and  simulated  several  composition  of  two  types  of 
active  materials  (LiNio.80Coo.15Alo.05O2  and  LiMn204).  They  pre¬ 
sented  that  their  blended  electrode  showed  combined  character¬ 
istics  (such  as  power  and  energy)  of  each  pure  active  material.  Jung 
and  Kang  [18]  introduced  a  multi-dimensional  model  of  lithium- 
ion  cells  with  blended-electrodes.  But,  their  model  was  an  empir¬ 
ical  model  which  could  not  be  used  to  predict  characteristics  of  a 
lithium-ion  cell  without  actual  experiment  data. 


This  paper  introduces  how  to  predict  performance  of  lithium- 
ion  batteries  with  blended-electrodes  by  mathematical  method. 
First,  we  introduce  the  blended-electrode  model  in  equilibrium 
state,  which  is  useful  to  predict  open-circuit  voltage  of  a  blended- 
electrode  with  arbitrary  combinations  of  multiple  active  mate¬ 
rials.  Second,  the  physics-based  dynamic  model  is  explained,  which 
can  simulate  battery  performance  in  dynamic  operating  conditions 
and  explain  physico-chemical  phenomena  inside  blended- 
electrodes. 

2.  Equilibrium  model 

A  pure  active  material  has  a  distinct  equilibrium  voltage  profile 
according  to  lithium  intercalation  level  or  lithium  stoichiometry, 
which  can  be  used  for  us  to  predict  charging/discharging  charac¬ 
teristics  of  a  cell  made  of  the  active  material.  By  combining  such 
equilibrium  voltage  data  of  several  active  materials,  it  is  possible  to 
predict  basic  characteristics  such  as  discharging  profile  of  a 
blended-electrode  as  follows. 

Fig.  1  represents  equilibrium  voltage  profiles  of  basic  active 
materials  for  this  study.  The  basic  concept  of  the  present  equilib¬ 
rium  model  of  a  blended-electrode  is  that  the  blended-electrode  is 
in  thermodynamically  equilibrium  state.  In  other  words,  the  surface 
energy  level  (or  the  potential)  of  the  constituting  active  particles 
should  remain  same  level  by  minimizing  potential  difference  be¬ 
tween  active  particles.  This  causes  each  different  active  particle 
may  have  different  state  of  charge  (lithium  stoichiometry)  under 
equilibrium  condition  of  the  blended-electrode.  Therefore,  we  can 
find  out  state  of  charge  (lithium  stoichiometry)  of  each  constituting 
active  material.  Theoretical  capacity  of  LiMn204  is  about 
150  mAh  g-1  while  usable  capacity  is  about  107  mAh  g-1  (150- 
107  =  43  mAh  g-1  is  reserved  to  protect  the  crystal  structure). 
Therefore,  minimum  lithium  stoichiometry  (soc  100%)  of  LiMn204 
is  stoich\oo  =  1.0  -  107/150  =  0.287  and  maximum  lithium  stoi¬ 
chiometry  (soc  0%)  is  stoicho  =  1.0  when  irreversible  loss  is  ignored. 
So,  we  can  setup  following  functions  using  Fig.  1  with  above 
relation. 


h(capai) 

f2  (stoichi)  where 
h(soci) 


stoichi  = 


capa[ 
capamax  i 


SOC; 


stoichi  -  stoich0  i 
stoich-ioo  i  -  stoicho  [ 


(1) 


Inversely,  lithium  stoichiometry  or  SOC  can  be  found  for  the 
given  equilibrium  potential  also  by  transforming  the  above  func¬ 
tions  as  follows: 


7,  '(U,)  =  cctpcii 
f2^(Ui)  =  stoichi  = 


capa  i 
capamax  i 


fpm 

capamax  i 


.  SOCj  - 

J  stoich^ooj  —  stoicho  i 


(ui)~  stoicho  j 
stoich'iooj  -  stoicho j 


(2) 


Note  that  the  maximum  capacity( capamax),  the  minimum 
lithium-stoichiometry  (stoicho),  and  the  maximum  lithium-stoi- 
chiometry(stoic/iioo)  of  the  active  material  i  are  known  values.  In 
order  to  get  a  full  voltage  profile  of  a  blended-electrode  with  n- 
constituting  active  materials,  we  calculate  accumulated  capacity, 
lithium  stoichiometry,  and  SOC  of  each  material  while  sweeping 
equilibrium  potential  from  llmax  to  Hmin  as  follows: 
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Capacity  [mAh/g] 
(a)  LMO  (LiMn2Q4) 
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Fig.  1.  Equilibrium  voltage  profile  of  pure 
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materials  used  for  the  present  study. 


capable nd  =  E  capaj 

1  =  1 
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stoichi  =  stoich[(U) 
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By  applying  the  scheme  explained  above,  equilibrium  potentials 
of  blended-electrodes  can  be  generated  as  shown  in  Fig.  2.  Positive 
blended-electrode  made  of  LiMn204  (briefly  ‘LMO’  from  now)  and 
Li(Nii/3Mni/3Coi/3)02  (briefly  ‘NMC’  from  now)  gets  more  capacity 
with  lowered  overall  potential  when  the  fraction  of  NMC  increases. 
Negative  blended-electrodes  made  of  graphite  and  soft  carbon  is 
simulated  as  shown  in  Fig.  2(b).  Soft  carbon  has  smoothly 
decreasing  profile  whereas  graphite  has  step-decreasing  profile 
due  to  lithium  staging  inside  crystal  layers  (see  Fig.  1).  Therefore, 
when  they  are  blended  together,  mixed  potential  profile  is  gener¬ 
ated.  As  the  present  blended-electrode  model  has  information  of 
each  constituting  material,  it  is  possible  to  extract  several  useful 


information  of  each  material  member  from  the  model.  Fig.  3(a) 
shows  the  contribution  of  each  material  member  to  the  whole  ca¬ 
pacity  production  by  the  blended-electrode  made  of  NMC  (70%) 
and  LMO  (30%).  Total  capacity  of  the  blended-electrode  is  about 
143  mAh  g-1.  NMC  and  LMO  produces  capacity  together  until  the 
net  capacity  reaches  65  mAh  g_1.  Flowever,  after  the  net  capacity 
reaches  65  mAh  g-1,  LMO  does  not  make  any  capacity  longer  while 
NMC  continues  to  produce  more  capacity  until  the  net  capacity 
reaches  143  mAh  g-1.  Fig.  3(b)  shows  different  view  of  capacity 
production  of  this  blended-electrode.  While  the  blended-electrode 
discharges  from  SOC  100%-0%,  SOC  of  LMO  component  already 
reaches  0%  when  SOC  of  the  blended-electrode  is  around  45%. 
Instead,  NMC  supports  capacity  production  to  the  last  of  the  dis¬ 
charging  process.  Note  the  SOC  slope  of  NMC  becomes  stiffer  after 
SOC  of  LMO  reaches  0%  in  Fig.  3(b)  since  NMC  solely  covers  the 
whole  capacity  production  from  this  point  on.  Lithium  stoichiom¬ 
etry  status  of  each  active  material  member  can  be  analyzed  as 
shown  in  Fig.  3(c).  Starting  from  0.3  at  SOC  100%,  lithium  stoichi¬ 
ometry  of  LMO  becomes  almost  1.0  when  SOC  of  the  blended- 
electrode  reaches  45%.  In  the  case  of  NMC,  lithium  stoichiometry 
changes  from  0.45  to  1.0  when  SOC  of  the  blended-electrode  moves 
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Fig.  2.  Simulated  equilibrium  potentials  of  blended-electrodes  according  to  blending 
ratio:  (a)  positive  electrode  composed  of  NMC  and  LMO,  (b)  negative  electrodes  made 
of  graphite  and  soft  carbon. 


from  100%  to  0%.  dQ/dV  profile  is  widely  used  for  analyzing  capacity 
production  of  an  electrode  according  to  voltage  window.  With  the 
present  blended-electrode  model,  one  can  separate  dQ/dV  of  each 
material  member  from  the  total  dQ/dV  profile  as  shown  in  Fig.  3(d). 
It  is  easily  identified  that  LMO  produces  capacity  in  the  voltage 
range  of  3. 8-4.2  V  whereas  NMC  covers  the  whole  voltage  window. 


solved  for  each  individual  active  particle.  In  addition,  it  is  consid¬ 
ered  that  those  particles  interact  with  each  other  to  simulate 
physical  behaviors  of  a  blended-electrode  as  shown  in  Fig.  4(b). 


3.1.  Formulation 

As  the  present  model  follows  the  framework  of  Newman  group’s 
[12,13],  this  paper  briefly  reviews  the  governing  equations.  Fig.  4(a) 
shows  that  there  are  three  mathematical  domains  in  the  model.  The 
following  electrolyte-phase  species  equation  is  solved  in  the  whole 
domain  (anode,  cathode,  and  separator). 


dfoe)  =  d_  f  neff^e\ 

at  ox  V  e  ax ) 


1  *+ju  where  D|ff  =  cTDe 
F  c 


(4) 


The  left  hand  side  term  is  the  transient  term,  and  the  right  hand 
side  term  consists  of  the  diffusion  term  and  the  source  term.  Here, 
Dff  represents  the  effective  lithium  ion  diffusivity  of  the 
electrolyte. 

The  electrolyte-phase  potential  equation  is  solved  in  the  entire 
domain.  Eq.  (5)  consists  of  the  ohmic  term  and  the  ion-diffusion 
term  with  a  source  term  representing  the  reaction  current  den¬ 
sity/1.  As  the  time  scale  of  the  potential  is  much  smaller  than  that 
of  Li-ion  species  equation,  the  transient  term  is  neglected  in  Eq.  (5). 


d_ 

Ox 


HD 


_0_ 

dx 


=  ~J 


;Li 


(5) 


It  is  known  that  the  ionic  conductivity  Keff  is  a  strong  function  of 
the  salt  concentration,  temperature,  and  material  composition. 
Following  relationship  is  used  for  finding  the  ionic  conductivity  of 
electrolyte  composed  of  EC  and  EMC  [19], 


k  =  -3.37115  +  12.5608m  -  7.89593m2  +  3.51922m3 

-  1.15471m4  +  18.1863X  -  6.22756 mx  -  13.6916m2x 

+  8.43904m3x  -  7.83732x2  +  19.607mx2  -  18.4529m2x2 

-  30.6369X3  +  29.2 mx3  -  0.0429918T  +  0.180877mT 

-  0.0836202m2T  +  0.0230098m3T  +  0.1959467x 
+  0.0676686 mTx  -  0.14134m27x  +  0.1474297x2 
+  0.173059m7x2  -  0.516347x3  -  0.000223097T2 
+  0.0001 1 1233m72  +  0.0000495286m272 

+  0.000952777T2x  +  0.001 1 7334 mT2x 

-  0.00061 91  57T2x2  -  3.46897  x  102773  -  2.75041 
x  1026 mT3  -  5.57653  x  102673x 

(6) 


3.  Dynamic  model 

Although  the  equilibrium  model  introduced  above  is  useful  for 
quickly  predicting  voltage  profile  and  analyzing  capacity  produc¬ 
tion  of  a  blended-electrode,  it  cannot  simulate  physical  behavior  of 
a  battery  under  dynamics  operating  conditions.  By  combining  the 
porous  electrode  theory  [12-14]  with  the  present  blended- 
electrode  model,  the  present  dynamic  blended-electrode  model  is 
developed.  Fig.  4(a)  shows  the  computation  domain  of  the  present 
model.  Finite  volume  mesh  is  constructed  in  x-direction  including 
negative  electrode,  separator,  and  positive  electrode.  Each  finite 
volume  has  multiple  active  particles  constituting  a  blended- 
electrode  system.  Spherical  lithium  intercalation  equation  is 


Note  that  m  is  electrolyte  concentration  and  x  is  fraction  of  EC  in 
the  above  equation. 

The  effective  ionic  conductivity  is  calculated  by  considering  the 
porosity  of  the  electrode  as  follows: 


Keff  = 


K 


eff 

D 


STK 

is  the  diffusive  conductivity  defined  as  follows  [20]: 


(7) 


,,eff 


2  RTKeff 


(8) 


The  solid-phase  potential  equation  (or  electron  transport 
equation)  which  follows  Ohm’s  law  is  solved  in  both  electrodes. 
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(a)  Contribution  of  NMC  and  LMO  to  capacity  production:  LMO  stops 
producing  capacity  when  the  net  capacity  production  of  the  blended-electrode 
reaches  65mAh  g'1  whereas  NMC  continues  to  produce  capacity  until 
the  net  capacity  reaches  the  maximum  value 


(b)  SOC  variation  of  each  active  material  member  (NMC  and  LMO) 
according  to  SOC  of  the  blended-  electrode:  When  the  net  SOC  reaches  40%, 
SOC  of  LMO  reaches  0%  whereas  remaining  SOC  of  NMC  is  still  60%. 


(c)  Li-stoichiometry  variation  of  each  active  material  member  (NMC 
and  LMO)  in  the  blended-electrode 


(d)  Capacity  production  curves  of  the  blended-electrode  and  the 
component  materials 


Fig.  3.  Simulated  result  of  a  blended-electrode  made  of  NMC  (70%)  and  LMO  (30%). 


ax 


=  r 


(9) 


where  creff  is  the  effective  electronic  conductivity  defined  as 
(7eff  =  (1  -  efa.  Again,  the  transient  term  is  ignored  in  the  above 
equation  where  the  electron  transport  is  very  fast  compared  to  the 
Li-ion  species  transport. 

As  explained  in  the  equilibrium  model,  open  circuit  potential 
(l/i)  of  each  active  material  is  defined  as  a  function  of  lithium 
stoichiometry  as  follows: 


Stoichj  =  (cs/cSim  ax)i  (10) 

Maximum  lithium  concentration  cs,max  is  theoretically  deter¬ 
mined  and  lithium  concentration  on  an  active  particle  cs  is 
calculated  by  solving  a  solid-phase  lithium  diffusion  equation.  In 
the  present  study,  all  active  particles  are  assumed  to  be  spherical. 


Therefore,  solid-phase  lithium  concentration  in  the  active  parti¬ 
cle  is  acquired  by  solving  a  spherical  diffusion  equation  as 
follows: 


^cs,i  _  j-v  92cs4  2Ds  i  0cs  i 
~dF  ~  s4  dr2  ~r  OF 


(11) 


The  above  spherical  equation  is  solved  for  each  particle  member 
constituting  the  blended-electrode.  Solid-phase  lithium  diffusivity 
in  active  particles  Ds  is  a  function  of  temperature  in  this  study.  Eq. 
(12)  shows  lithium  intercalation  speed  will  increase  at  high  tem¬ 
perature,  which  may  result  in  improved  C-rate  performance  of  a 
cell. 


Dsj  =  Dr5fexpl^f 


(12) 
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Fig.  4.  Schematics  of  the  present  dynamic  blended-electrode  model:  (a)  Mathematical  domain  of  the  quasi-2D  blended-electrode  model,  (b)  Inter-particular  relationship  among 
three  distinctly  different  active  particles  located  in  a  single  mathematical  mesh  domain. 


Loading  fraction  ^  of  each  active  material  in  a  blended- 
electrode  is  expressed  as 

'/'i  =  771  where  wa  =  V«a,i  (13) 

Wa  itl 

Average  density  of  active  materials  is 


exp 


^a,i^  „ 

~w 


Vi  -exp  - 


~RT 


Vi 


where  V\  =  4>s  ~  <t>e  ~ 


(17) 


Exchange  current  density  i0,i  which  is  a  function  of  lithium 
concentration  in  electrolyte  phase  and  solid  phase  is  expressed  as 
follows  [14]: 


■  1st 


(14) 


Therefore,  porosity  of  an  electrode  which  has  active  material, 
binder,  and  conductive  agent  can  be  calculated  as  follows: 


-  1  -  Ge)  -  Ge)  -  Ge)  (,5) 

Specific  interfacial  area  (as)  which  reflects  the  effective  reaction 
area  of  each  active  material  is  defined  as  follows: 


A- 

asi  &  Nu1-  where 
L± 


'  A[  =  47 zrf 
Vi  =  1 7 w? 

N-  -^L 

1  ftVi 


(16) 


Reaction  kinetics  at  the  interface  between  each  active  particle 
surface  and  electrolyte  is  assumed  to  follow  the  Butler— Volmer 
equation  as  follows: 


*o,i  —  [&(Ce)aa  (Cs,max  _  cs)  a(cs)Q!cji.  (^8) 

When  a  blended-electrode  consists  of  two  different  active  par¬ 
ticles  (binary  system),  two  different  active  particles  will  face  each 
other  and  interactive  reaction  will  occur  to  build  equilibrium  be¬ 
tween  them.  The  driving  force  of  this  interactive  reaction  is  voltage 
difference  between  two  active  particles.  If  active  particle  1  (i  =  1) 
and  active  particle  2  (z  =  2)  constitute  a  blended-electrode,  the 
interactive  reaction  current  density  can  be  defined  as  follows: 

Jl2  =  fcl2(as,12'o,12)AL/l2exp|^^T--^| 

where  AUU  =  U,  -  U2  (19) 

In  Eq.  (19),  as<u  is  specific  surface  area  between  active  particles  1 
and  2,  which  is  assumed  to  be  the  smaller  value  between  two 
particles.  In  the  same  way,  exchange  current  density  between  two 
particles  (i0, 12)  is  determined,  ku  is  reaction  coefficient  between 
two  particles,  which  is  an  assumed  value  in  the  present  study  since 
it  is  difficult  to  measure  this  value.  For  each  active  particle, 
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vLi  _  vLi  I  vLi 
J 1  —  Jrxn,l^Jl2 
vLi  _  vLi  I  vLi 
J 2  ~  J rxn,2  +  J21 


(20) 


In  each  electrode,  total  reaction  current  is  the  sum  of  reaction 
current  from  each  active  material  member.  Therefore,  the  following 
description  can  be  made. 


ju  =  =#+#  (21) 

i  =  1 

This  total  transfer  current  density  is  plugged  into  the  governing 
equation  Eqs.  (4),  (5),  (9)  and  (11). 

When  the  blended-electrode  is  a  ternary  system  which  has 
three  distinctly  different  active  particles,  interactive  reaction  cur¬ 
rent  densities  among  particles  are 


In  the  same  way,  the  present  model  can  be  extended  to  any 
blended  electrode  system  with  N-different  active  particles. 

As  the  physical  behavior  of  a  cell  is  strongly  affected  by  tem¬ 
perature,  it  is  necessary  to  find  the  cell  temperature,  especially 
under  high  C-rate  case  where  it  is  difficult  to  maintain  constant 
temperature.  It  is  assumed  that  temperature  in  the  cell  is  same 
everywhere  due  to  the  small  scale  of  the  present  model.  Therefore, 
the  following  lumped  heat  equation  is  coupled  with  the  present 
model. 


mcP  fa  —  hAsmf(T  T'amb)  +  Q  (24) 

Heat  source  term  q  consists  of  the  followings  [21]. 


J12  =  ^12(as,12^o,12)^^12-exp|^^  } 

J23  =  ^23  (as,23  *0,23)  At/23  'exp{“f^  ^7^  _  fj  | 
J31  =  fc3l(as,3l'o,3l)At/3iexp{^^T--lj| 


where  A Uu  =  -  U2 

where  AL/23  =  U2  -  U3 
where  A l/31  =  U3  - 


(22) 


<7+  =  E  (asio)i(0s  -  -  Ui)  +  E  (Osio)i  (T§).+  ffeffV2<As  +  KeffV20e  +  Vln  CeV0e 

i  i  v  7 1 

qs  =  KeffV20e  +  KQfVln  CeV^e  (25) 

q~  =  E  (Osio)i (<^s  —  —  Uj)  +  E  (Usio)i(r .  +  creffv20s  +  Keffv2(fie  +  Kpffvln  cev</)e 

i  i  v  7 1 


Current  density  of  the  blended-electrode  becomes 

[  JV  =  Jrin.l  +Jl2  +Jl3 
jU  =  Y) J'P  =  Jl1  +J2  +J3  where  {  H  =  Jrxn,2  +^21  +J23 

1  =  1  [M  =  J™,3 +J31  +J32 

(23) 


Note  each  electrode  has  reaction  heat,  electronic  ohmic  heat, 
ionic  ohmic  heat  whereas  separator  has  only  ionic  ohmic  heat  in 
the  above  expression. 

Governing  equations  with  boundary  conditions  for  the  present 
model  are  summarized  in  Table  1.  The  model  is  incorporated  into 
FORTRAN  and  solved  by  finite  volume  method. 


Table  1 

Governing  equations  and  boundary  conditions  of  the  Li-ion  cell  model. 


Description 

Conservation  equation 

Boundary  condition 

Charge  (electrolyte  phase) 

Charge  (solid  phase  in  electrode  matrix) 

A  (*'*$)+£(*?£>»*)  =-JLi 
£(«**&)  =JLi 

90e  _  0<^>e  1  _  f) 

dx  X=Q  dx  \X=L 

_n-eff0^s  1  _  ffCfityA  _  L 

dx  x=0  +  dx  \x=L  A 

Species  (electrolyte  phase) 

Species  (solid  phase  in  particles) 

9cs  rj  92cs  ,  2 Ds  9cs 

0t  —  Us  dr2  _l_  r  0r 

rt-effa&l  _  neffd<t>s\  _  n 

-  9x  \x_Ln  ~  °+  dx  x_in|s  —  U 

0Ce|  _  0Ce|  _ 

dx  x=0  —  0x  lx=I 

cs  =  cs  o  at  t  =  0 

Ds^~  =  0  at  r  =  0 
or 

|Li 

Ds— -  =  J—  on  particle  surface 
or  r 

Energy 

mcPW  =  hAsmf(T  -  Tamb)  +  <j 
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Table  2 

Simulation  result  of  NMC/LMO  blended-electrode  according  to  the  material 
composition. 


Total  loading 
(mg  cm'2) 

NMC  loading 
(mg  cm'2) 

LMO  loading 
(mg  cm'2) 

NMC 

fraction  (%) 

Electrode 
capacity 
(mAh  cm  2) 

16.0 

0.0 

16.0 

0 

1.630 

16.0 

4.8 

11.2 

30 

1.863 

16.0 

8.0 

8.0 

50 

2.096 

16.0 

11.2 

4.8 

70 

2.212 

16.0 

16.0 

0.0 

100 

2.445 

(a) 


Fig.  5.  Simulated  C-rate  performance  of  NMC/LMO  blended-electrode:  (a)  Tce n  =  298  K, 
(b)  TCeii  =  273  K. 


Table  3 

Active  material  composition  of  full-cells  constructed  for  the  present  study. 


Full-cell  no. 

Cathode 

Anode 

NMC  (%) 

LMO  (%) 

Graphite  (%) 

Carbon  (%) 

1 

70.0 

30.0 

90.0 

10.0 

2 

50.0 

50.0 

90.0 

10.0 

3 

30.0 

70.0 

5.0 

95.0 

3.2.  Result  and  discussion 

LMO  and  NMC  are  widely  used  active  materials  for  Li-ion  cells. 
LMO  has  advantage  of  fast  lithium  intercalation  due  to  spinel 
structure  along  with  relatively  high  operational  voltage,  which 
enables  a  Li-ion  cell  to  produce  high  power.  For  this  reason,  LMO 
will  be  a  good  choice  as  positive  material  for  HEV  application.  In 
contrast,  NMC  has  advantage  of  larger  theoretical  capacity 
(277  mAh  g-1)  compared  to  LMO  (150  mAh  g_1)  although  lithium 
intercalation  rate  is  a  bit  lower  compared  to  LMO.  Therefore,  it  is 
reasonable  to  think  that  a  blended-electrode  composed  of  both 
NMC  and  LMO  has  intermediate  characteristics  between  them.  In 
other  words,  it  is  possible  to  acquire  both  advantages  from  NMC 
and  LMO  with  a  blended-electrode. 

Discharging  performance  of  NMC/LMO  blended-electrodes  is 
simulated  using  the  dynamic  model  according  to  five  blending 
cases  as  listed  in  Table  2.  Calculation  result  shows  net  electrode 
capacity  linearly  increases  according  to  increment  of  NMC  fraction. 
This  result  is  almost  identical  with  the  result  of  equilibrium  model 
(see  Fig.  2(a)).  Effect  of  blending  ratio  on  C-rate  capability  is  shown 
in  Fig.  5.  It  is  clear  that  blending  more  LMO  leads  to  improved  C-rate 
capability  due  to  fast  lithium  diffusion  in  LMO  particles,  especially 
under  high  C-rate  condition  over  7C.  In  addition,  LMO  helps  a  cell 
produce  more  power  or  large  current  at  low  temperature. 

In  order  to  validate  the  present  dynamic  model,  three  different 
full  cells  with  different  electrode  compositions  are  built  as  listed  in 
Table  3.  Baseline  model  parameters  used  for  the  simulation  are 
listed  in  Table  4.  Three  cells  are  simulated  under  the  given  test 
conditions,  and  compared  with  experimental  result  as  shown  in 
Fig.  6.  By  blending  different  active  materials,  each  cell  shows 
distinctly  different  voltage  profiles.  Cell  #1  and  cell  #2  have  most  of 
their  capacity  production  in  the  voltage  window  of  4.1 -3.4  V 
whereas  cell  #3  shows  slowly  decreasing  voltage  profile  between 
4.1  and  2.8  V.  This  result  came  from  anode  material  selection 
(graphite  vs.  amorphous  carbon).  As  cell  #1  and  cell  #2  share  the 
same  anode,  difference  between  them  comes  from  the  cathode 


Table  4 

Baseline  model  parameters. 


Parameter 

Unit 

Value 

Note 

Radius  of  NMC 

pm 

4.5 

Measured 

Radius  of  LMO 

pm 

7.5 

Measured 

Radius  of  graphite 

pm 

7.5 

Measured 

Radius  of  soft  carbon 

pm 

3.8 

Measured 

Li-ion  diffusivity  of  NMC 

cm2  s'1 

5.5  x  10'10 

Fitted 

Li-ion  diffusivity  of  LMO 

cm2  s'1 

2.5  x  10'9 

Fitted 

Li-ion  diffusivity  of  graphite 

cm2  s'1 

4.5  x  10'10 

Fitted 

Li-ion  diffusivity  of  soft  carbon 

cm2  s'1 

1.7  x  10'9 

Fitted 

Thickness  of  negative  current 

pm 

10.0 

Measured 

collector 

Thickness  of  positive  current 

pm 

20.0 

Measured 

collector 

Thickness  of  negative  electrode 

pm 

50.0 

Measured 

Thickness  of  positive  electrode 

pm 

55.0 

Measured 

Thickness  of  separator 

pm 

25.0 

Measured 

Porosity  of  negative  electrode 

0.28 

Measured 

Porosity  of  positive  electrode 

0.25 

Measured 

Porosity  of  separator  electrode 

0.37 

Measured 

Ionic  conductivity  in  electrolyte 
Transference  number 

S  cm'1 

Eq.  (6) 

0.363 

From  ref.  [6] 

Faraday’s  constant  (F) 

C  mol'1 

96,487 

Li-ion  diffusivity  in  electrolyte 

cm2  s'1 

2.5  x  10'6 

From  ref.  [6] 

Average  concentration  of  Li-ion 

M 

1.0 

in  electrolyte 

Matrix  electronic  conductivity 

S  cm'1 

0.062 

Measured 

(positive) 

Matrix  electronic  conductivity 

S  cm  1 

0.1 

Measured 

(negative) 

Votage  IV]  Voltage  [V]  Voltage  IV] 
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0  500  1000  1500  2000  2500  3000  3500  4000 

Time  [sec] 

(a)  Discharging  profiles  of  full-cell  #1 


0  1000  2000  3000  4000  5000  6000  7000  8000 

Time  [sec] 


(b)  Discharging  profiles  of  full-cell  #2 


Time  [sec] 

(c)  Discharging  profiles  of  full-cell  #3 
Fig.  6.  Experimental  validation  of  Li-ion  cells  with  blend  ed-electrodes  (full-cell  #1,  #2,  #3 ). 


side.  Discharging  simulation  of  IOC  was  successful  with  cell  #2 
which  has  large  LMO  fraction  whereas  cell  #1  showed  maximum 
7C  discharging  performance  with  the  present  model.  This  result 
confirms  that  using  much  amount  of  fast  diffusion  material  in  the 
electrode  is  beneficial  for  getting  high  power  from  a  lithium-ion 
cell  although  there  are  other  factors  to  control  C-rate  perfor¬ 
mance  of  a  cell  such  as  electrode  thickness.  Cell  #3  achieved  very 
high  C-rate  performance  (20C)  compared  to  other  cells.  This  high 
power  performance  was  possible  with  the  following  three  reasons. 
First,  the  cathode  consists  of  large  amount  of  manganese  spinel  and 
the  most  of  active  material  fillings  the  anode  is  amorphous  carbon 
which  has  fast  solid-phase  lithium  diffusion  speed.  Second,  elec¬ 
trode  thickness  of  cell  #3  is  less  thick  (about  60%)  compared  to  cell 
#1  and  #2.  Third,  both  electrode  of  cell  #3  is  very  porous  compared 
to  the  other  cells,  which  means  amount  of  active  material  loading  is 
low.  Consequently,  energy  density  of  cell  #3  is  lower  than  the  other 
cells.  In  contrast,  opposite  approaches  of  the  above  three  methods 
should  be  considered  when  a  cell  with  high  energy  density  is 
designed.  High  energy  density  cells  can  be  applied  to  energy  stor¬ 
age  systems  (ESS)  which  does  not  require  high  power  capability  as 
xEVs  do.  For  automotive  applications,  cell  #1  and  #2  will  be 


Normalized  length  [x/L] 

(a)  Negative  (Graphite)  ||  Positive  (NMC) 


Normalized  length  [x/L] 

(b)  Negative  (Soft  carbon)  ||  Positive  (LMO) 

Fig.  7.  Li-stoichiometry  profiles  of  each  active  material  in  full-cell  #1  during  3C 
discharging. 
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appropriate  for  PHEV  or  EV,  whereas  cell  #3  will  be  a  good  choice 
for  HEV  which  require  high  power  than  large  amount  of  energy. 

It  is  possible  to  analyze  physical  behavior  of  each  active  material 
member  in  a  cell  with  the  present  model.  Fig.  7  show  lithium 
stoichiometry  of  each  active  material  in  cell  #1  under  discharging 
condition.  It  is  observed  that  lithium  stoichiometry  window  of 
graphite  is  wider  than  that  of  carbon  in  the  anode  while  lithium 
stoichiometry  window  of  LMO  is  wider  than  that  of  NMC  in  the 
cathode. 

When  a  blended-electrode  has  several  constituting  active  ma¬ 
terials  with  quite  different  physical  properties  such  as  temperature 
dependency  of  lithium-ion  diffusivity,  the  voltage  profile  may 
significantly  change  according  to  the  operating  condition.  Fig.  8 
presents  simulated  discharging  profiles  of  a  lithium-ion  cell  with 
a  blended-positive  electrode  (NMC  85%  and  LPF 15%)  and  a  negative 
electrode  (graphite  100%).  Under  the  mild  operating  condition 
(room  temperature,  298  I<),  discharging  ends  up  with  negative- 
limiting  condition  and  the  positive  voltage  plateau  of  LPF  shows 


Fig.  8.  Simulated  1C  discharging  profiles  of  a  Li-ion  cells  with  positive  electrode 
composed  of  NMC  85%//LFP  15%  and  negative  electrode  composed  of  graphite  100%,  (a) 
T  =  298  K,  (b)  T  =  273  K. 


up  lately  (see  Fig.  8(a)).  On  the  other  hand,  the  positive  voltage 
plateau  of  LPF  shows  up  early  and  the  discharging  ends  up  with 
positive-limiting  condition  when  the  cell  operates  at  low  temper¬ 
ature  (see  Fig.  8(b)).  In  this  simulation  case,  it  was  assumed  that  the 
temperature  dependency  of  lithium-ion  diffusivity  in  NMC  particle 
is  stronger  than  that  of  graphite  particle.  Therefore,  solid-phase 
mass-transport  in  NMC  particle  is  getting  slower  than  that  of 
graphite  particle  at  low  temperature.  As  the  lithium  concentration 
on  the  surface  of  NMC  particle  is  quickly  saturated  due  to  slow 
lithium  intercalation,  operation  of  LPF  starts  up  earlier  under  the 
low  temperature  condition  than  the  mild  operating  condition.  In 
other  words,  the  length  of  span  A  and  B  strongly  depends  on  the 
material  properties  of  the  constituting  active  particles,  which  can 
be  easily  predicted  with  the  present  dynamic  model. 

Another  simulation  examples  of  lithium-ion  cells  with  blended- 
electrodes  are  shown  in  Fig.  9.  According  to  the  blending  ratio  and 
negative/positive  ratio,  numerous  combinations  of  cells  can  be 
constructed.  When  material  properties  of  each  active  material  are 
prepared,  characteristics  and  physical  behavior  of  those  cells  can  be 


(a)  Discharging  profiles  of  a  Li-ion  cell  with  positive  electrode  composed  of 
NMC  50%  //  LMO  50%  and  negative  electrode  composed  of  graphite  90%  // 
LTO  10% 


(b)  Discharging  profiles  of  a  Li-ion  cell  with  positive  electrode  composed  of 
NMC  70%  //  LFP  30%  and  negative  electrode  composed  of  graphite  100% 

Fig.  9.  Simulated  result  of  various  Li-ion  cells  with  blended-electrodes. 


194 


S.  Jung  /  Journal  of  Power  Sources  264  (2014)  184-194 


predicted  with  the  present  model.  This  will  eventually  expedite  cell 
development  process  which  usually  takes  much  time  and  effort. 

4.  Conclusion 

Two  mathematical  models  of  lithium-ion  cells  with  blended- 
electrodes  are  developed.  First,  the  equilibrium  model  can  simu¬ 
late  any  combination  of  active  materials  in  a  blended-electrode  by  a 
quick  and  easy  method.  The  equilibrium  model  enables  us  to 
decompose  a  blended-electrode  into  each  active  material, 
analyzing  steady-state  characteristics  of  each  constituting  member 
such  as  voltage  window,  and  lithium  stoichiometry.  Second,  the 
physics-based  dynamic  model  can  conduct  simulation  under 
various  operating  conditions  such  as  C-rate  and  temperature  by 
solving  physico-chemical  governing  equations.  Simulation  results 
show  good  agreement  with  the  experiment  data,  confirming  that 
the  present  model  is  useful  for  checking  candidate  active  material 
compositions  at  the  early  stage  of  battery  development. 

Nomenclature 


Acronyms  and  abbreviations 

as  active  surface  area  per  electrode  unit  volume  (cm-1) 

A  area  (cm2) 

c  concentration  of  lithium  (mol  cm-3) 

cp  specific  heat  (J  K-1) 

capa  capacity 

D  diffusivity  of  lithium  (cm2  s-1) 

Fact  activation  energy 

h  heat  transfer  coefficient  (W  cm-2  K-1) 

z'o  exchange  current  density  (A  cm-2) 

I  applied  current  (A) 

j  volumetric  current  density  (A  cm-3) 

k  reaction  coefficient 

L  thickness  (cm) 

m  mass  (g) 

N  number  of  active  particles 

r  radial  coordinate,  radius  (cm) 

R  universal  gas  constant  (8.314  J  mol-1  K-1) 

soc  state  of  charge 

stoich  lithium  stoichiometry 

t  time  (sec) 

t°+  transference  number  of  lithium  ion  with  respect  to  the 

velocity  of  solvent 
T  temperature  (I<) 

V  volume  (cm3) 

x  thickness  coordinate  (cm) 

U  open  circuit  potential  (V) 

q  heat  source  (W) 

Greek  letters 

aa,  olc  anodic  and  cathodic  transfer  coefficient  for  an  electrode 
reaction 
e  porosity 


0  potential  (V) 

7]  surface  overpotential  of  an  electrode  reaction  (V) 

k  ionic  conductivity  (S  cm-1) 

kd  diffusive  conductivity  of  a  species  (A  cm-1) 
p  density 

a  electronic  conductivity  (S  cm-1) 

t  Bruggeman  tortuosity  exponent 

0  loading  fraction 

w  loading  amount  (g  cm-2) 

?  capacity  fraction 

Superscripts  and  subscripts 
a  active  material 

b  binder 

amb  ambient 

c  conductive  agent 

diff  diffusion 

e  electrolyte  phase 

eff  effective 

eq  equilibrium 

i  ith  member 

Li  lithium 

ref  reference 

rxn  reaction 

s  solid  phase 

surf  surface 
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